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1. Introduction. 



The dynamical symmetry breaking (DSB) has been considered as a pos- 
sible mechanism of appropriable fermion mass generation in the quantum 
field theory for a long time [1- 20]. Usually, the broken symmetry is the 
chiral one, which plays a very important role in the high energy physics [3- 
12, 18, 19]. There are some methods of DSB investigation: the direct cal- 
culation of composite fields effective action [1, 6, 7, 16], Schwinger- Dyson 
equations (SDE) analysis [3- 8, 13, 15, 17, 18] or solution of composite field 
renormalizations group equations [20] . 

These methods have been generalized for curved spacetime case in 4D 
QED [21], four- fermionic models [22- 28], Einstein general relativity [29] 
and 2D induced gravity [30, 31]. The results of the DSB in curved space- 
time have been the subject of considerable attention because the curvature 
induced phase transitions, accompanied by the creation of non- zero vac- 
uum expectation values of both elementary [32, 33] and composite bosonic 
[33- 35] or fermionic [21- 33] fields, is turned out to be important to con- 
struct the realistic Early Universe scenario. 

It is well-known that, unfortunately, the self-consistent quantum theory 
of gravity does not exist now. Because the Einstein general relativity is 
unrenormalizable [36- 37] , we have to find out the other models with more 
attractive quantum features. The most natural and simplest generalization 
of Einstein theory is the R 2 - gravity [33, 38- 44]Q. It is renormalizable and 
asymptotically free theory. However, it has some disadvantages, such as the 
non-unitarity in the usual pertubative theory sense and unphysical ghosts 
presence. From the modern point of view the higher- derivative gravity is 
nothing but the model including the next term in the low enrgy expansion 
of future complete quantum theory, based perhaps, on (super)string theory. 
This opinion lets us don't take into account the last problems because we 
anyway have to work below the Planck scale, naturally limiting the low 
energy gravity physics [33] . 

Therefore, the investigation of the DSB both in 2D induced in the (su- 
per) string theory R 2 -gravity [45- 48] and in the 4D original version of this 
model provides us some essential information about the possible features 
of future complete quantum theory of gravity. 

1 For the review of recent results, see [33]. 
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In the present paper we apply SDE formalism to quantum R 2 gravity 
with the fermions on the flat background. The ladder or rainbow approx- 
imation, when the vertex of fermion- graviton interaction and graviton 
propagator are taken to be free is used. Both Landau- like general covari- 
ant and conformal gauges are considered in 2D case. The covariant gauge 
providing the minimal structure of graviton Green function (GF) is chosen 
for 4D spacetime. SDE are obtained and the integral equations determing 
the exact fermion GF are written evidently. The numerical analysis of 
their solutions is done in details. The DSB is shown to exist in a differ- 
ent gauges. The dependance of dynamical fermion mass on the coupling 
constant is found out. 



2. Dynamical symmetry breaking in two dimensional 
R 2 - gravity. 

We will consider here the theory with the following action: 

1 

L2M 2 



S = I d l x^~g —(R 2 + 2A) + i^(x)D^ , (1) 



where R is the space-time curvature, ip- 2D spinors, A- cosmological 
constant and the dimensional value M is the fermion- graviton coupling 
constant. 

D, = 0, + l -ufa ah (2) 

is the spinor covariant derivative with the standard spin- connection u^. 

Local Dirac matrices 7 M (#) can be expressed through the usual flat ones 
7 a and tetrads e%(x) : 7 M (x) = e^{x)^ a and, finally, a a b = |[7a,7&]- Greek 
and Latin indices correspond to the curved and flat tangent spacetimes 
respectively. 

In accordace with the general background field method one gets: 

9iiv = + h^y. (3) 
We choose the gauge- fixing term in the folowing form [49] : 

S gf = 7=§ ^ 2 x v ^(V A ^-/? 2 V^)(^V p V^ 3 V^V^)(V (T ^-/3 2 V^), 

(4) 
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where h = and fli, (32,03 ~ gauge- fixing dimensionless parameters. The 
most useful tool for the calculation of the graviton Green functions for 
these kinds of complicate expressions is the projective operators method 
[50]. It gives 



} 1 k " ( j]^k v k G + rf a k v k p + rfWk* + rf'WlA + 



A(A - A/c 4 ) 

+ ^(A-l)(ft-l/2)(l + A) _ _2_) + 

2 + a _ mzjrzim yr _ * + ^)] , (5) 



where 

S = -P1W2 ~ 1) 2 (1 + Ps)k 8 + -f- (l + 4ft (ft - l/2) 2 (/? 3 + 1)) • (6) 
Graviton - fermion interaction vertex has the usual form [29, 30]: 

IVO, fc) = ^(2p A + k x )Y(^r]x<jr]^ - m^V™ ~ m^^)- (7) 

The most general Lorentz invariant form for the inverse exact fermion 
GF is the following: 

S-\p) = A(p 2 )tf-B(p 2 ), (8) 

where A(p 2 ) and B(p 2 ) are the unknown functions we should find out. 
The SDE for GF (8) in the ladder approximation are given by: 



S- 1 -S Q 1 (p)= [ AQr^(q,p-q)S(q)r Xa (p,q-p)G^(p-q) (9), 

where So(p) = 1/ p 1 is the free fermion GF. 

These equations determine the extremum of composite fermionic fields 
effective potential [16]: 

V ef f = -iSp [ ln{S Q l S) - S^S + 1 ] + V 2 , (10) 
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Vo = 



where 

1 f d 2 p 



d 2 q 

2 J J2nf] (2^ 



Sp [T^p - q, q)S(q)r pa (q-p 1 p)S(q)\ G^(p - q) 

(11) 

corresponds to the two-particle irreducable vacuum diagrams. Formula 
(11) is written in the ladder approximation, where the vertices and graviton 
GF are chosen to be free and the only fermion GF is taken to be exact. 

Now we have all the necessary parts of Feynman diagrams to calculate 
the final expressions for the SDE (9) and the effective potential (10). After 
the Wick rotation and tedious algebra the integral equations for structure 
functions A{x) and B{x) are obtained: 



A(x) = 1 + 9 ( 
Jo 



A(y)dy 



1 



K A {x,y) 



r yA\y)+B\y)x 

B{x) = g L STM %,y) 

and the effective potential (10) is given by 



(12), 
(13) 



T 2 ,1 

V '" = ~ii\Jo 



ln(A 2 (x) + 
V x 



B 2 {x)\ A{x)(a{x)-1)+B\x) 
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dx 



dy 



l A 2 (x)x + B 2 (x)J A 2 {y)y + B 2 {y) 

B(x)B(y)K B (x,y) 
where T is the ultraviolet cut off parameter, 



A 2 (x)x + B 2 (x) 
A(x)A(y)K A (x,y) + 



+ 



(14) 



1 M 2 P 2 Q 2 a i 2a a, ^ x B (P 2 ) 

9 = X = T*' V = h A{P ) = A{X) ' B{X) = T ■ ( 5) 

The evident form for the integral equations kernels KA(x,y), KB(x,y) 
have been obtained explicitly for the arbitrary (3i, fy, 03- However they 
are too large to present them here. 

We would note only that these kernels contain, in general, the terms 
with unpleasant factors like (x — y)~ l . It means that the infrared diver- 
gencies caused by the virtual massless gravitons whose momentum tends 
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to zero appear here. This type of divergences arises , probably, because 
the quantum corrections for the interaction vertex are omitted in our ap- 
proximation [21, 30]. 

However in the conformal and Landau- like gauges, which have the 
most considerable physical meaning, these divergences don't take place 
[30, 31]. That is why we investigate here the f3\ — > oo case, corresponding 
to Landau- like gauge. Then 



K A {x,y) = 



ft xy + {x + y){x + y + h)l± (2ft - l/2)(x - y)' 



(ft " I) 2 



{x + y + l 2 )-^-j^{{x-yf-ll) 



yj(x + y + l 2 ) 2 - 4xy 



+ 



Pi xy+{x + y){x + y- / 2 )/4 (2/? 2 - 1/2) (x - y)' 



(ft - I) 2 
(x - y) 2 



K B (x,y) = 



h /i(ft-l/2) 
&(x,yjl /2 ) , 2fa{x + y) 



(x + y-l 2 ) + 



+ 



ft 2 



2(ft - l) 2 



y/(x + y-l 2 ) 2 -4xy /i(ft - 1/2) 
x + y\ , {x-yf i B(x,y, 4 /2 ) 



a: (16) 



1-2 



ft 2 



L2(/? 2 - l)^ 



„x + y\ (x — y) 
1 + 2—^ + 



/ 5 



/ 2 ; ' h(p 2 -l/2)\^/( x + y -l 2 )2-4xy 
2 ^ 1 



h(P 2 - 1/2) J ^(x + i/ + / 2 ) 2 - 4x1/ 
— 2|x — z/| 



h<Jh-l/2)' 



(17) 



where /i = ^, / 2 = v^i 



/3 2 -l/2 



and 



e(x, 2/, «) = e( (Vi - Vv) 2 - « 2 ) - 0(a 2 - (v^ + W) 2 )- (is) 

The analytical solution of thes non-linear integral equations does not 
seem to be possible. Therefore, we present here the results of numerical 
calculations by means of the standard iterative procedure, described, for 
example, in [29- 31] 

The dependence of structural functions A{x) and B(x) on the Euclidian 
momentum square is presented in the Fig. 1 for the different values of 
coupling constant g and fixed l\ = 4 and (3 2 = 1.05. A{x) (curve 1) doesn't 



almost depend on g. The only trivial solution B = exists for function 
B(x) for small g. However for g > g c = 0.23 the type of solution changes 
essentially and only the non-trivilal ones increasing with g growth provide 
the minimum of effective potential (14). The curves 4, 3, 2 correspond to 
the following values of g = 0.25, 0.30, 0.35 respectively. 

The dependence of dynamical mass, defined by the pole of exact fermion 
GF, on the gauge coupling constant g is shown in the Fig. 2. It means 
that after the analytical continuation into the pseudoeuclidian region the 
dynamical mass is obtained as the solution of the equation 

m 2 A(m 2 ) - B(m 2 ) = 0. (19) 

This Figure shows us clearly the typical behaviour of order parameter m 2 
in the course of phase transition accompanied by the creation of bifermionic 
condensate in the above critical region g > g c . 

Let us discuss now the conformal gauge, which plays the very important 
role in the (super)string theory. In this gauge 

g^y = exp((p)r]^. (20) 
Then, the graviton propagator is given by: 

M 2 

G(k) = ^ (21) 
and the fermion- graviton interaction vertex are the following [31]: 

r(p,fc) = i(2|f+ (22) 

The equations for the structural functions and and effective potential 
have the same form as in the previous case (12)- (14) with: 

„ ( N x 2 + y 2 + 6xy - l{x + y) Q(x,y,l 1 / 2 ) 
K A {x,y) = 



4/ y/(x + y - I) 2 - 4xy 

x 2 + y 2 + 6xy + l(x + y) 1 

4/ y/(x + y + l) 2 -±xy 

2x + 2y-l e(x,y,lV 2 ) 
K B (x,y) = ==+ 

21 v( x + y - l ) 4x v 



(23) 
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2x + 2y + l 1 

2/ ^/{x + y + l) 2 -4xy 

where / = — ^ > 0. 

The plot of the functions A{x) and B(x) are presented in the Fig. 3. 
The dynamical symmetry breaking takes place in this case as well. The 
behaviour of function A(x) doesn't almost depend on the value g (curve 1). 
The non-trivial solutions for function B(x) appear for g > g c = 2.0 only. 
The curves 4, 3, 2 correspond to the g = 3,4.5,5.5. Figure 4 represents 
the dynamical mass as a function of coupling constant g. 

3. Schwinger- Dyson equations in 4D higher- deriva- 
tive gravity. 

The similar program can be realized for the four- dimensional gravity with 
the Lagrangian, containing square curvature terms ( for details, see book 
[33]): 



S g = J (fx^g~ [ai R 2 + a 2 R^R^ + a z R llV p a R' wpa + Aj , (25) 

directly following [33, 49, 50]. ai,«2,«3 are the arbitrary constants here. 
We choose the parameters of gauge fixing action: 

S » / = § / d " a; ^ r 9( V ^i-AV ( ,ft)W"'V,V' > +/3 3 VV)(V„^-/3 2 V^), 

(26) 

demanding that the square on h^ v part of action, (S + S g f)^ won't contain 
the non-minimal terms. Then: 

{S g + S gf )W = I f tfxy/^h^frforrripan 2 ~ \Pi&%U 2 + 



2 „ 

2 V 2 

where 



+ «(«fTri P * - 6%) W\ (27) 



&% = \W + S W) (28) 
and 

R a n 4«i+a 2 2ai + a 2 + 2a 3 

p 1 = a 2 + 4a 3 , fo = 77 r, = — ; (29) 

4(ai - 03) 02 + 4^3 



Then, the graviton GF is given by: 
1 



Qpvpo 



A(/c 4 + 2A)L 



2((3 2 k 4 + AW 
/c 4 (l - 4&) - 2A 



(30) 



where A = The graviton- fermion interaction vertex has the same 
form as for the covariant gauge in the 2D model (7) and we obtain the 
same expressions for the SDE (12)- (13) and, with the accuracy up to the 
positive constant factor, for the effective potentilal (14) but, of course, 
different kernels: 

1 



Ka(x, y) 



y 



+ 



X 



2Zi 



+ y + k ~ y/(x + y + k) 2 -Axy) 
+ J/)) + 



16 32x 

l{ - xy + 7(x + yf/i 11 



+(x + y-k - \/(x + y- k) 2 - 4xyG(x, y, l\ /2 )) x 
xy + 7/4(x + y) 2 11 



2k 



--g-0* + J/))+ 
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4(1 - 4/? 2 ) + 27 + /2 ~ V / (^ + ?/ + ^) 2 -4x ? /) x 

~^ + J/)) + 



^2xy + l/2(x + ?/) 2 1 



/ 5 
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4(1 - 4&) 



(x + y-h - \/(x + y-l 2 ) 2 -4xyG(x, y, l\' 2 )) x 



2.-»/ + l/2(x + ;/) 2 1 

+ 2 ( * + 2/) 

9 r / 

K B {x,y) = U8l 2 x [k ((2x + 2y - k) (x + y - k - 



(31). 



128/fx 

^{x + y-kf-Axy^x, y, l\' 2 ))- 
(2x + 2y + k)(x + y + k ~ ^ {x + y + k) 2 - Axy^j - 

l 2 [{2x + 2y-l 2 ){x + y-l 2 - ^(x + y - l 2 ) 2 - 4xyS(x, y, / 2 1/2 ))- 

{2x + 2y + l 2 ){x + y + l 2 - ^{x + y + / 2 ) 2 - 4^))] , (32) 
where / 2 = , q = ^ > 1/4> p = (le^)- 1 . 
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On the Fig. 5 the plot of functions A(x) (curve 1) and non-trivial 
solutions of B(x) minimized the effective potential for g > g c = 2.5 (curves 
2, 3, 4 corresponds to g = 4.5,4.0,3.4) are presented. = 2, l\ = 0.1 
here. Fig. 6 represents the dymamical mass dependence on the coupling 
constant g. 

4. Conclusions 

Thus, the quantum higher-derivative gravity interacted with the fermions 
has been shown by the numerical analysis to contain the phase with the 
dynamically broken chiral symmetry and dynamically generated fermionic 
mass. The most important result of our paper is the fact that existance 
of such phenomena doesn't depend on the gauge choice in the different 
models of R 2 - gravity. However the values of g c and the character of phase 
transition are strongly depend on both the parameters of the model and 
the gauge type. 
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